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Abstract. We investigate a ring of N FitzHugh-Nagumo elements coupled in phase- repulsive 
fashion and submitted to a (subthreshold) common oscillatory signal and independent Gaussian 
white noises. This system can be regarded as a reduced version of the one studied in [Phys. Rev. 
E 64, 041912 (2001)], although externally forced and submitted to noise. The noise-sustained 
synchronization of the system with the external signal is characterized. 
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INTRODUCTION 

In Ref. 03 — through comparison with the synchronization patterns arising in two- 
dimensional arrays of FitzHugh-Nagumo (FHN) elements with phase-repulsive linear 
nearest-neighbor coupling — the authors were able to conclude that intracellular cal- 
cium oscillations in cultures of human epileptic astrocytes do interact, since the phases 
of nearby oscillating astrocytes maintain a nontrivial relationship. It is a fortunate fact 
that the (space-independent) FHN model is one of the very few multicomponent systems 
for which a nonequilibrium potential (NEP) has been found |2l|3], since NEPs allow in 
general for a deep insight on the dynamical mechanisms leading to pattern formation and 
other phenomena where fluctuations play a constructive role [4J. The (albeit minimal) 
extension of the result in Refs. [0 [3j towards extended systems carried out in this work 
is however enough to shed light on the dynamical cause of the conclusion in Ref. 0]|: a 
dynamical symmetry breakdown takes place because the phase-repulsive coupling min- 
imizes the corresponding NEP. When the system is externally forced with a frequency 
less than the typical inverse deterministic time the cycle duplicates, breaking down into 
an "excited" phase and an "inhibited" one. These phases force neighbor elements to al- 
ternate with the one in between, thus creating a nontrivial phase relationship between 
nearby oscillating elements. 

The system we consider is sketched in Fig. [1] a ring of N = 256 identical FHN el- 
ements with phase-repulsive nearest-neighbor coupling and submitted to a (subthresh- 
old) common oscillatory signal and independent Gaussian white noises £ m (t), <3v, W 
with (% m (t)% n (t')) = 2r\8 mn 8{t — t'), m,n = l,...,2iV. The set of equations governing 
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FIGURE 1. Sketch of the system and of its response A c {t). 



its dynamics is 

iii = a c Ui(l-uf)-Vi + S g (t)-D(ui + i+Ui-i) + ri^ Ui (t) + r2^ Vi {t) (1) 
Vi = e(^Ui-Vi + C) + r 3 ^ u ,(t)+r4^ Vi (t), i=l,...,N, u N+ \=u\. 

where £ = 0.01 is the ratio between the relaxation rates of w, and v,-, /3 = 0.01, a c = 0.06 
and C = 0.02 is a suitable constant to set the rest point in Fig. [2^. D = 0.01 is the 
phase-repulsive coupling constant, and the r ; (which determine the transport matrix) are 
n = 0.998 x 10 2 , r 2 = 0.499 x 10 1 , r 3 = 0.998, r 4 = 0.499 x 10" 1 . Moreover, taking the 
Milshtein integration step as dt = 5 x 10 -3 , we estimate the typical inverse deterministic 
time as 0.838 x 10~ 3 and so we take the excitation frequency Q.q as a fraction of that 
value (typically 0.1-0.4). Given that, S g (t) = 0.0275sinQ ^ 

THE NONEQUILIBRIUM POTENTIAL 

Excitable dynamics can be conceptually decomposed into two phases, a fluctuation- 
dominated one and a deterministic one. It would be highly desirable to find a Lyapunov 
function, since it greatly simplifies the dynamical analysis. However, the existence of 
non- variational (or conserving) components in the phase- space flow is a hint that the 
integrability conditions fail for the purely deterministic system. This apparently insur- 
mountable drawback was partially solved two decades ago by Graham and collaborators 
(see references in [j21) who defined the NEP for Langevin-type dynamics as the zero- 
noise limit of the logarithm of the stationary probability density function (pdf). The 
extra freedom in the choice of the transport matrix can render in some cases the prob- 
lem integrable. That is precisely the case for the space-independent FitzHugh-Nagumo 
model in its bistable and excitable regimes OEl. This approach can be generalized to 
extended systems and the NEP associated to Eq.Q (in the adiabatic limit, i.e. for slow 



FIGURE 2. a) Phase-space excursions in excitable regime (the nullclines are indicated in dashed line); 
b) Time evolution of u (dashed line), v (dotted line), and the NEP (full line) during a phase-space 
excursion. The scales of v and the NEP were adjusted for better comparison. 




FIGURE 3. a) NEP in excitable regime; b) Stationary pdf in excitable regime. 
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FIGURE 4. a) NEP in bistable regime; b) Stationary pdf in bistable regime. 
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FIGURE 5. Time evolution of u for two neighbor neurons. 

which must obey the integrability condition /3Ai + fa/s = 2X [O. 

Figure [2j) depicts (in full line) the time evolution of the NEP during the phase- space 
excursion starting at the upper initial condition in Fig.Ek together with that of u (dashed 
line) and v (dotted line). We remark that in Figs. [2^ and[2j) there is no noise and 4>(f) is 
the Lyapunov functional of the deterministic dynamics. Figures [3^l and [3]) (respectively 
[3k and[4j)) are 3D and contour plots of the NEP and the corresponding stationary pdf for 
the excitable (respectively bistable) regime. 

RESULTS FOR THE COUPLED SYSTEM 

Synchronization between the coupled system and the external signal is observed above 
some noise-intensity threshold. Figure [5] is a plot of the time evolution of m (full line), 
together with that of iij+i (dashed line) for a given neuron i, showing their phase relation 
to the signal (dotted line^j According to Fig. [2j we may call "active" those cells i for 
which Ui(t) exceeds some threshold value u±. Because of the coupling, as one neuron 
becomes active, it inhibits the activation of its nearest neighbors. The perfect alternance 
seen in the figure may fail because of the noise, a necessary ingredient for the activation. 

A detail of the alternance can be seen in Fig.|6]for an = 21 subset of the ring. Figure 
[6k shows a situation (snapshot) of poor synchronization, in which only two neurons are 
active; Fig. [6}) exhibits a case of a "kink" in the synchronized configuration, induced by 
the fact that noises are local. Note that the kinks break locally the observed coherence, 
and the complete history of the time evolution can be followed as a record of activity 
(see Fig. [7j>. 



The signal has been augmented in about two orders of magnitude and shifted to aid the sight. 
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FIGURE 6. Two snapshots of {m,} showing different degrees of synchrony. 




A measure of "activity" for the whole ring is 

1 N 

Ac W = 77£0k-O)-"th]- (3) 
iV (=i 

In perfect synchrony, Ac = 0.5. Note that since the signal is subthreshold for the coupled 
system, Ac = below threshold. Figure [8^ depicts the activity as a function of time 
for a fixed noise intensity, showing again its phase relationship with the signal (dashed 
line). In Fig. [8j) we show the NEP for the whole ring as a function of time, together 
with the (scaled) signal for reference. We remark that the observed dynamical symmetry 
breakdown decreases the Lyapunov function of the whole ring with respect to that of the 
homogeneous state, providing the route to stable synchronization. 
A global estimator of synchronization can be defined as 

fj Ac(t)dt 

G » = \^r- (4) 

Figure [9^ is a plot of G a as a function of the noise intensity. The existence of a 
threshold value of noise intensity and of a saturation effect can be clearly seen. The 
noise intensities are low enough not to degrade the excitable dynamics. 
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FIGURE 8. a) Ac vs / for high synchronization; b) Time evolution of the NEP 





FIGURE 9. a)G a vsTj;b)Cvsrj. 



Numerical simulations indicate that the coherence of firing decreases with the noise 
intensity although the global activity (representative of global estimators) keeps the 
order of magnitude. To quantify this phenomena we have calculated the normalized self- 
correlation C = (uiUi + 2) as a function of the noise intensity ?). As we show in Fig. [9}) 
the system shows a kind of "stochastic resonance in coherence" that cannot be inferred 
from measures of global activity. 



CONCLUSIONS 



We have investigated the noise-induced synchronization with an external signal of a ring 
of phase-repulsively coupled FHN elements. We have derived the exact NEP of the ex- 
tended system and the observed symmetry breakdown was related with the Lyapunov- 
functional properties of the NEP. We remark that the same conclusion holds qualita- 
tively for the work in Ref. 0]. Although the observed phenomenon is noise-sustained 
and global activity increases with noise intensity, a degradation of coherence can be 
appreciated. 
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